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Abstract. Consider a discrete uniformly elliptic divergence form equation on 
the d dimensional lattice Z'* with random coefficients. In [3] rate of conver- 
gence results in homogenization and estimates on the difference between the 
averaged Green's function and the homogenized Green's function for random 
environments which satisfy a Poincare inequality were obtained. Here these re- 
sults arc extended to certain environments with long range correlations. These 
environments are simply related via a convolution to environments which do 
satisfy a Poincare inequality. 



1. Introduction. 

In this paper we continue the study of solutions to divergence form ehiptic equa- 
tions with random coefficients begun in [3] . In [3] we were concerned with solutions 
u^x^rjjio) to the equation 

(1.1) riu{x,ri,uj) + V*a{Tri:Uj)Vu{x,ri,u}) = h{x), x G Z*^, w <S fi, 

where Z'' is the d dimensional integer lattice and {^l,J-',P) is a probability space 
equipped with measure preserving translation operators t^; : f2 — > f2, x G Z"^. The 
function a : ft ^d{d+i)/2 ^^.^^-^ q ^j^^, 

space of symmetric dxd matrices satisfies 

the quadratic form inequality 

(1.2) Xld < a(cj) < A/rf, uj en, 

where Id is the identity matrix in d dimensions and A, A are positive constants. 

It is well known [71 [TOl [H] that if the translation operators t^, x S Z'^, are 
ergodic on H. then solutions to the random equation (|1.1[) converge to solutions of a 
constant coefficient equation under suitable scaling. Thus suppose / : R"^ — > R is a 
C°° function with compact support and for e satisfying < e < 1 let Me (a;, 77, w) be 
the solution to ()1.1|1 with h{x) = e^f{ex), x £ Z''. Then Ue{x/e,e'^r],uj) converges 
with probability 1 as e — > to a function Uhom(2;, rj), x G R'^, which is the solution 
to the constant coefficient elliptic PDE 

(1.3) ?7Whom(2:, Jy) - VahomVuhom(a::, 77) = /(x), x € R'', 

where the dxd symmetric matrix ahom satisfies the quadratic form inequality 
(|1.2p . This homogenization result can be viewed as a kind of central limit theorem, 
and our goal in [5] was to show that the result can be strengthened for certain 
probability spaces {fl,J-,P). In particular, wc extended a result of Yurinskii [T2] 
which gives a rate of convergence in homogenization: 
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Theorem 1.1. Let f : R'^ — > R &e a C°° function of compact support, Ue{x,ri,ijj) 
the corresponding solution to U.l\) with h{x) = e'^f{ex), x E Z'', and Uhom(a;, x G 
R'', the solution to il.S]) . Then for certain strong mixing environments {^l,J-,P)- 
see discussion below- there is a constant a > depending only on d, A/A and a 
constant C independent of e such that 

(1.4) sup ( |we(a;/£,£^7?, •)-ithom(a;,77)l^ ) < C'e", for < e < 1. 

It is evident that any environment $7 can be considered to be a set of fields 
w : Z'^ — > R" with n < d{d + l)/2, where the translation operators r^, x £ Z'', 
act as TxU!{z) = uj{x + z), z G Z'', and a(tj) = a(a;(0)) for some function a : 
R" R'i('i+i)/2. Yurinskii's assumption on (17, J", P) is a quantitative strong 
mixing condition. Thus let x{-) be a positive decreasing function on R+ such that 
limg^.oo x('z) — 0- The quantitative strong mixing condition is given in terms of 
the function x(-) as follows: For any subsets A,B of Z*^ and events Ta, Tb C ^, 
which depend respectively only on variables i-u(x), x € A, and w(y), y B, then 

(1.5) |p(r^nrB)-P(r^)P(rB)| <x( inf 

In the proof of (|1.4p he requires the function x{') to have power law decay i.e. 
limg_j.oo q°'x{l) = for some a > 0. Evidently (jl.Sp trivially holds if the w(x), x g 
Z'^, are independent variables. Recently Caffarelli and Souganidis [2] have ob- 
tained rates of convergence results in homogenization of fully nonlinear PDE under 
the quantitative strong mixing condition (jl.Sp . In their case the function xil) 
is assumed to decay logarithmically in q to 0, and correspondingly the rate of 
convergence in homogenization that is obtained is also logarithmic in e. In their 
methodology a stronger assumption on the function x(-), for example power law 
decay, does not yield a stronger rate of convergence in homogenization. 

In [9] Naddaf and Spencer obtained rate of convergence results for homoge- 
nization under a different quantitative strong mixing assumption than (|1.5|) . Their 
assumption is that a Poincare inequality holds for the random environment. Specif- 
ically, consider the measure space (17, J^) of vector fields a; : Z'' — )■ R'^, where F 
is the minimal Borel algebra such that each uj{x) : O — >■ R'^ is Borel measurable, 
X G Z'*. A probability measure P on {(l,J^) satisfies a Poincare inequality if there 
is a constant m > such that 

(1.6) Var[G(-)] < X;{\\dcjG{-)f ) for aU functions G : f7 ^ C, 

where da,G{y,uj) = dG{ui)/doj{y), y £ Z"*, is the gradient of G(-). In [£! it is 
assumed that P is translation invariant i.e. the translation operators Tx, x £ Z**, 
acting by Txijj{z) = uj{x + z), z G Z'^, are measure preserving, and that the Poincare 
inequality (|1.6p holds. Rate of convergence results are then obtained provided 
sJyuj) = a(w(0)) in (dH), where the function a : R*^ ^- R'*('i+i)/2 is G^ and has 
bounded derivative, in addition to satisfying p.2p . In recent work Gloria and 
Otto ini [5] have developed much further the methodology of Naddaf and Spencer, 
obtaining optimal rates of convergence in homogenization under the assumption 
that the environment satisfies a weak Poincare inequality. This weak Poincare 
inequality holds for an environment in which the variables a(r2-ti;), x £ Tfi , are 
independent, whereas the inequality ()1.6|) in general does not. 
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In this paper we shall obtain rate of convergence and related results for homog- 
enization in environments defined by a(aj) = a(w(0)) where uj : Tl'- ^ R" is a 
convolution u:{-) = di S fi. The function h : TJ^ ^ R" ® R*^ from TJ^ to 

n X k matrices is assumed to be q integrable for some q < 2, and the probability 
space {^l,J-, P) to satisfy the Poincare inequality ()1.6|) . Unlike environments which 
satisfy a Poincare inequality, these environments a;(-) can easily be shown to have 
long range correlations. In particular, consider the massive field theory environ- 
ment studied in [3] consisting of fields : Z'' — > R with measure P formally given 

by 



(1.7) exp 



- 1/ (V0(a;)) + im20(a;)2 



JJ^ (i(/i(a;)/normalization. 



where V : R^ — > R is a uniformly convex function. Then {Vt,J-,P) with measure 
(jl.7|) satisfies the inequality p.6p . In the Gaussian case when V{-) is quadratic one 
has that the correlation function ( (f){x)(t){Q) ) = Gm2(a;), x £ Z'^, where the Green's 
function Gn{-) is the solution to 

(1.8) r)G,-i{x)+V*V"VG^{x) ^ 5{x), .t e Z"*. 

Hence the correlation function ( 0(x)0(O) ) decays exponentially in |x| as |x| — > oo. 
Taking a;(-) — h*- </>(•) for some h S L''{'L'^) we have that 

(1.9) ( w(a;)w(0) ) - ^ /i(a; - y)h{^y')G,r.2 {y - y') , 

and so if 1 <q< 2 then {uj{fff) < oo. Choosing now ^.(z) = l/[l+\z\'^/'^+^], z e 
for any e > we see from (|1.9p that ( uj{x)uj{0) ) ~ jxl^^"^ as |x| — > oo. 

The limit as m — >■ of the measure (|1.7p is a probability measure P on gradient 
fields w : Z'' — )- R'', where formally oj{x) = V(/)(a;), a; € Z'^, a result first shown 
by Funaki and Spohn [4|. This massless field theory measure satisfies a Poincare 
inequality (|1.6p for all d > 1. In the case d = 1 the measure has a simple structure 
since then the variables cj(a;), a; S Z, are i.i.d. For d > 3 the gradient field theory 
measure induces a measure on fields : Z'^ — > R which is simply the limit of the 
measures (|1.7|) as to — >■ 0. For d = 1,2 the m ^ limit of the measures (|1.7I) on 
fields (/> : Z^ ^ R does not exist. If d > 3 then ( (t){x)4:{0) ) ~ \x\-^'^^'^'> as oo 
for the massless field theory. Observe now that 

(1.10) 0(x) = [VGo(.T-y)]^V(^(j/) = /i*w(x), xeZ'', 

where Go(-) is the Green's function for with -q = 0,V" = Id- Since h : Z'' ^ K'' 
in (jl.lOp is g integrable for any q > d/{d — 1), the environment of massless fields 
: Z'' — )■ R with d> 3 satisfies the condition g < 2 of the previous paragraph. 

Rather than attempt to formulate a general theorem for environments uj = h*uj 
where {(l,T,P) satisfies the Poincare inequality (|1.6p . we shall only rigorously 
prove that the results obtained in [3] hold for massless fields i/i : Z'^ — )• R with 
d > 3. In §2 we indicate the generality of our argument by showing that the proof 
of Proposition 5.3 of [3| formally extends to environments ui = h * Cj. In §3 we 
implement the method of §2 to prove the following theorem for massless fields: 
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Theorem 1.2. Let a : R — > ^d{d+i)/2 ^ (ji j^jiqUqji qji yj^fi values in the 
space of symmetric d x d matrices, which satisfies the quadratic form inequality 
and has hounded first derivative Da{-) so ||_Da(-)||oo < oo. For d > 3 let 
{fl,J-,P) be the probability space of massless fields determined by the limit 
of the uniformly convex measures Jj. 7| ) as m — > 0, and set a(-) in to be 

a((/)) = a(0(O)), G 51. Then Theorem 1.1 holds for the probability space {fl,J-,P). 
Let Ga_^jj{x), X e Z'*, be the averaged Green's function for the random equation lll.l\) 
and Gai„„,77(x), x £ Hf^, the Green's function for the homogenized equation il.3\} . 
Then there are constants a,7 > depending only on d and the ratio A/ A of the 
constants A, A of hl.2jl . and a constant C depending only on ||Z)a(-)||oo, A/A, d such 
that 

(1.11) |Ga.,(x) - Ga,„„.,(x)| < -^^^-^^^^e-'v^l-l X e Z'^ - {0}, 

(1.12) |VG,.,(x) - VG,,_.,(x)| < -^^^^-^p-^e'^V^I-l, x^Z'^- {0}, 

(1.13) |VVGa.,(x)-VVGa,„.„,,(:r)| < A(|:,| + l)^+^ "''^'^' ^■eZ'^-{0}, 
provided Q < i] < K. 



2. Variance Estimate on the Solution to a PDE on n 
Following §5 of [3] we consider the solution $(^, rj^uj) to the equation 

(2.1) 77$(^,,,,c^) + a|a(c^)a5<I>(e,?7,^) = -a|a(c^), 77 > 0, ^ G ^ e fl, 

and let P denote the projection orthogonal to the constant function. Then our 
generalization of Proposition 5.3 of [3] is as follows: 

Proposition 2.1. Suppose a(-) in \2.1\) is given by a{uj) = a(cL>(0)) where a : R" — 
j^d{d+i)/2 ig a d X d symmetric matrix valued function satisfying the quadratic 
form inequality lll.2\) and ||Da(-)||oo < 00. The random field u) : Z"^ ^ R" is a 
convolution uj{-) = h*u}{-) of an nxk matrix valued function /i : Z"^ — > R'^^R*^ and 
a random field uj :T/' ^ R'^ . The function h is assumed to be po summable for some 
po with 1 < Po < 2 and the probability space {fl,J-,P) of the fields a) : Z'' — > R'^ 
to satisfy the Poincare inequality \1.6\) . Then there exists pi depending only on 
d,A/X,pQ and satisfying 1 < pi < 2, such that for g e LP{Z'^,C'^ (g) C*) with 
1 < P < Pi and v £ C^, 

(2.2) \\P J2 g{x)d,Htv,r.-)v\\ < ||fe||poHgHp , 

where C depends only on d, n, fc, A/A,po- 

Proof. From (|1.6p we have that 
(2.3) 

IIP 9i^)dim,ri,r.-)vr < ^ W a|i) ^ 5W5e$(^, r?, r.-)«lP • 
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From the chain rule we see that 
(2.4) ' 



duj{y) 



h{y - z) 



Hence using the translation invariance of the probability measure P on we con- 
clude from (1231), (123) that 
(2.5) 



d 



''duj{y) 



We define a function u : Z"* x ^ C*^ by 

(2.6) u[z,u) = e-''< K'^{y,i,V,r,Lo)v]h{y + z) , 

yezd 

where d^^{-,£_,7],uj)v : Z^ — >• C" is the gradient of ^{S,,r],uj)v with respect to 
a; G il. Observe now that 

d 



h{y - z) 



(2.7) Vu(a;-z,a;) = e*("-^'-« ^ 



'aa;(2/) 



h{y - z) 



whence (|2.5p becomes 
(2.8) 

IIP ^ ff(x)9^$(^,,;,r,.)«|P < E II E .9(a:)e^(^-^>«Vzi(x - z, •) |p . 

a;GZ<i zeZ"* ajeZ^* 

Next we take the gradient of equation (|2.ip with respect to w(-). Using the 
notation of [3] we have that 

(2.9) ri d^<^{y, v, + Dla{uj{0))D^ d^<^{y, V, 

= -Dl[5{y)Da{u;iO)){v + d^^{tv,^)v}] foryGZ^wG^!. 



Evidently (|2.9p holds with w G replaced by t^uj for any z G Z''. We now multiply 
(j2.9p with r^cj in place of uj on the right by e~^^'^h{y + z) and sum with respect to 
y G Z''. It then follows from ([^ that 

(2.10) ?7 u(z,a;) + V*a(w(z))Vu(z,w) = -V*f{z,Lu), 
where the function / : Z'' x ® C'^ is given by the formula 

(2.11) f{z,uj) = Dk{io{z)){v + d^^{i,r^,T,Lo)v}e-''<h{z) . 

Now a5$(e, ?7, •)« G L2(n, C"^) and \\d^<^{S,, v, < A\v\/X. Hence iihe L'^{Z'^, R" 
R*^) then the function / is in L^iZ"^ x n^C^ (g) C'') and ||/||2 < || £>&(•) ||oo(l + 
A/A)|t;|||/i|j2. We see from fU^i that if / G ^^(Z^ x fi, ® C*^) then Vu is in 
L2(Z'^ xn^C^® C^) and || Vu||2 < H/lb/A. It follows then from and Young's 
inequality that (|2.2I) holds with po = 2 and p = 1 . 

To prove the inequality for some p > 1 we use Meyer's theorem [8]. Thus for 
any 1 < g < oo we consider the function / as a mapping / : Z'' L'^{fl, (g) C*"') 
with norm defined by 

(2.12) ||./||^ = J2 ll/(^'-)ll 
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where ||/(z, OH 2 is the norm of /(z, •) e L'^{n,C'^ (g)C''). Noting that the Calderon- 
Zygmund theorem apphes to functions with range in a Hilbert space [11], we con- 
clude that there exists qo depending only on d, A/A with 1 < go < 2 such that if 
ll/llgo < 00 then II Vw||g < 2||/||g/A for qo < q < 2. If ft, is po integrable with po < 2 
we can take maxIpoi'Zo] = 9i < 9 < 2. It follows again from (|2.8p and Young's 
inequality that ()2.2p holds with pi = 2qi/{3qi — 2). □ 



3. Proof of Theorem 1.2 

It will be sufficient for us to establish the conclusion of Lemma 5.1 of [3] for the 
massless field theory environment {il,T,P) of Theorem 1.2. Using the notation of 
[3] we have the following: 

Lemma 3.1. Let {^1,J-',P) be an environment of massless fields : Z"^ — > R 
with d > 3, and a : R ^ ^d(d+i)/2 

as in the statement of Theorem 1.2. Set 
SL{(j)) = a((/)(0)), (/) € ri. Then there exists pq{A/ X) with 1 < po{A/X) < 2 depending 
only on d and A/ A, and a constant C depending only on d such that 

(3.1) WTr.Jp.o. < _ A/A)-/2 for 1 < p < po(A/A) . 

Proof. It will be sufficient for us to bound HT^.^gHoo in terms of H^Hp for g : Z'^ — 
(g) C'' of finite support. Let Q be a cube in Z'^ containing the support of the 
function g(-) and {Q,q,Tq, Pq^m) be the probability space of periodic functions 
(i : Q — > R with measure 



(3.2) exp 



-^F(V^(x)) + im20(x)2 



(i(/)(a;)/normalization, 



where we assume m > and V : R'' ^ R is with a(-) = V"{-) satisfying 
the quadratic form inequality (|1.2p . We denote by Qq the space of periodic fields 
(D : Q — )■ R'' and let F : J7q x fig — )• C be a function which for some constants 
A, B satisfies the inequality 
(3.3) 

\F{Cj,^)\ + \d^F{y,Co,4>)\^d^F{y,i::,ct>)\ < Aexp[B{||tD||2+||0||2}] , y & Q, Co ^ Uq, ^ ^ Vie, 

Letting (•)Q,m denote expectation with respect to the measure p.2p we see from 
the Brascamp-Lieb inequality [T] that the Poincare inequality 

(3.4) VarQ,™[^^(V0,0)] < \{ ||d^F(V0, </))|p )q.„ + A( \\d^F{V <i>,<j>)f )q.„ 

A 

holds. We shall show using p.4p that ||T,._^(7||oo is bounded in terms of \\g\\p if the 
environment is the probability space {flQ,J-Q,PQ jn)- The result will then follow 
by taking first Q Z'^ and then m 0. 

Let us suppose that the cube Q is centered at the origin in Z'^ with side of 
length L, where L is an even integer. Let G,, : Z'^ — > R be the solution to (|1.8p 
with V"{-) = Id and Gjj^q : Q R the corresponding Green's function for the 
periodic lattice Q, so 

(3.5) G^q(^) = ^G„(.T + in), xeQ. 
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Then any periodic function cf) : Q H can be written as 

(3.6) <Pix) ^ J2^VG^,Qix-y)fV^{y) + Y,7iGr,M^-y)(t>iy) , x e Q. 

yeQ y£Q 

Taking 77 = in (|3.6p we have a representation 

(3.7) (/){■) ^ hi + kg ^ <!){■), where ||/iq||,< for g>d/(d-l), 

and llfcgllg < C,/min[L''(i-i/«),i2] for g > 1 and q 7^ (i/(d - 2) . 

In (j3.7p the vector /ig = ['^Q,i) ^Q.d] is a column vector, the operation * denotes 
convolution on the periodic lattice Q, and Cq is a constant depending only on q, d. 

We first prove (j3.f p when r = 1. For the environment {ilq, J^q, PQ,rn) we have 
from dSIll) that 

(3.8) ri,,5(C,0) = ^.g(a;)e-"-«Pb(/j5*i;(a:) + fcQ*0(a:)) . 

Let 7im{Q) be the Hilbert space of functions / : Hq — > which are square 
integrable with respect to the measure Pq^m- It follows from p.4p that if u g C'' 
the norm of Ti rjgiS., ■)v S T~Lm{Q) is bounded as 
(3.9) 

m^r,9{t -H' < ^ E E II E g{x)hQ^,{^^z)Dh{m)vf+^ E n E 5(a;)fcQ(x-z)i?b(0(x))Hp . 

Since d > 3 we can choose q such that d/{d—l) < q < 2 and q^ d/{d — 2). It then 
follows from ^1}, that for p = 2q/{'iq - 2) > 1 

(3.10) ||Ti,,g(e,-)t'|P < CM\l\\DH-)\\Uv? 

where a{q) = 2min[(i(l — l/(7),2]. Let (il, J^, P™) be the probability space of fields 
: Z'^ R with measure Pm given by (|1.7p . Proposition 5.1 of [3] enables us to 
take the limit of p.lOp as Q Z*^ to obtain the inequality 

(3.11) ||ri,,.9(^,.)«lP < cM\l\\DH-)\\loHV>^ 

for the environment (fi, Pm). Finally Proposition 6.1 of [3] enables us to take 
the limit of (|3.1ip as m — >■ provided d > 3. We have proved p.ip when r = 1. 

To prove the result for r > 1 we consider the environment (QqjTq, PQ^m) and 
write as in [3] 

(3.12) Tr,,Mi,cj,{-))v = PE5(2^)e-'"-^b(0(x))95P,(r,0(.)) , 0(-) € 

xeQ 

where the functions Fr{4){-)) are defined inductively by 

(3.13) IPriH-)) + dld^FMi-)) = Pdl[h{m)diFr-i{H-))l r > 2, 

It is easy to see that d^Fr e 'Hm{Q) and 

(3.14) \\d(,Fr\\ < (1-A/A)'->| forr>2. 
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Using the representation p.7p for (/)(•) we can consider the Fr, r > 2, defined by 
p.l3p as functions of Cj{-) and (j>{-), which we denote by Fr{ui,(j)). Observe now 
that for I < j < d, 

(3.15) -^^5(2^)e-"-«b(/i§*(I;(.T)+fcQ*0(a;))a^F,(T,w,T,0) = 

xeQ xeQ 
Let j : Q x fig — > C be given by the formula 

(3.16) u.,,(z, </>(•)) = e-'^-^Y.^'t'Fr{y,r.cl}{-))hQ^,{y + z) . 
Then as in ^M, (|2Jl) we have that 



(3.17) Vw,,,(a;-z, </)(•)) = e^^^"-)'* J] r_,-|^a5F.(r,0(.));iQ,,(y - z) 

= e'(--)-«r_,— ^95F,(r,<:^,r,(/)) . 

()UJj{z) 

Similarly to (|2.10p . (|2.1ip we see that Urj(z,0(-)) satisfies the equation 

(3.18) ^ ?.,,,(z, </)(•)) + V*V^.,.,(z,^(-)) = FV*/r.,(2>(-)) , 
where the function f^.j : Q x VIq — > is given by the formula 

(3.19) /2,,(z» = Dh[c^{z))ve-''<hQ^,{z) , 

UA^A) = £'b((/)(^))a5F,_i(r,0(.))e-'^-«/iQj(z)+b(0(z))V?/,_ij(zX.)), r >2. 
From (|3TT5l) . (|3T7)) we sec that 

(3.20) -||-_-^5(x)e-'-«b(/i5*L^(x) + fcQ*0(x))a^F,(T,ci,T,0)f < 

II ^ g(a;)e-"-«/iQ,j(a;-z)Z?b(0(a;-z))a5^^,(T,_,0(.))f + 11 ^ 5(x)e-*^-«b(0(x-z))Vu^ 

Observe now from (|3.14p and Young's inequality for functions with values in a 
Hilbert space that 

(3.21) II E gi^y-'^'-^hQAx. - z)Dh{4>{x - z))d^Fr{T,^,,cj,{-))\\^ 

z^Q xeQ 

1 2 



< C 



\\Dhi-)\U\gUhQj,il-X/AY-'\v 



where p = 2q/{3q — 2) with 1 < q < 2 and C depends only on d. We can bound the 
second term on the RHS of p.20p similarly. Thus let L'^{Q,'Hm{Q)) be the Banach 
space of functions f : Q 'Hm{Q) with norm 

(3.22) 11/11^ = 5]||/(x)r . 

xeQ 

From ([339]) it follows that 
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(3.23) ||/2,,||, < Cpb(.)!|oo||/iQ,,i|,|i'| , 

WfrJU < q|i^b(.)||oo||/iQj||,(l-A/Ar-2|i;| + (l-A/A)||V^.,_i,,||, if r > 2, 

where C depends only on d. We see from the Hilbert space version of the Calderon- 
Zygmund theorem (see [TT] page 45) apphed to p.lSp that for q > 1 there is a 
constant (5(g) > such that 

(3.24) W^urjq < [l + 5iq)]\\frjq andhm%)=0. 
It foUows then from (jX^ . ([5:^ that 

(3.25) lUJ, < Cr\\Dh{-)\U\hQj,[l + S{q)r-^l-X/AY-'\v\, 
where C depends only on d. Now Young's inequality and (|3.24p . p.25p imply that 

(3.26) II g{x)e-''<h{<f>{x - z))Wurj{x - z, ?!)(-))f 

2 



< C 



rPb(.)||oo||.g||p||/»Q,,||,[l + s{q)r-\i ~ x/Ar-'\v\ 



where p = 2q/{3q — 2) with 1 < g < 2 and C depends only on d. We can argue 
now as in the r = 1 case to establish p.ip for r > 2 by choosing g < 2 to satisfy 
[1 + %)](!- A/A) < (l-A/A)i/2. □ 
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